Current-driven quantum criticality in itinerant electron ferromagnets 
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We determine the effect of an in-plane current flow on the critical properties of a 2d itinerant 
electron system near a ferromagnetic -paramagnetic quantum critical point. We study a model in 
which a nonequilibrium steady state is established as a result of exchange of particles and energy 
with an underlying substrate. The current j gives rise not only to an effective temperature equal to 
the voltage drop over a distance of order the mean free path, but also to symmetry-breaking terms 
of the form j ■ V in the effective action. The effect of the symmetry breaking on the fluctuational and 
critical properties is found to be small although (in agreement with previous results) if rotational 
degrees of freedom are important, the current can make the classically ordered state dynamically 
unstable. 

PACS numbers: 73.23.-b,05.30.-d,71.10.-w,71.38.-k 



One of the most intriguing issues in modern condensed 
matter physics is quantum criticality, the complex phe- 
nomena associated with phase changes driven by varia- 
tion of Hamiltonian parameters such as pressure, mag- 
netic field, or chemical composition in the limit as tem- 
perature T — > 0, 0, 0, B| ■ Quantum phase transitions 
differ from classical phase transitions because the non- 
commutativity of position and momentum in quantum 
mechanics implies that spatial and temporal fluctuations 
are coupled at a quantum critical point. While quantum 
phase transitions in equilibrium systems have been the 
subject of extensive study, phase transitions caused by 
nonequilibrium drives (for example an imposed current 
flow) have been much less extensively studied. Scaling 
theories and calculations have shown that one important 
effect of a departure from equilibrium is to produce de- 
coherence effects which may be modeled as an effective 
temperature @, @, 0, H, HJ- However, dep artures from 
equilibrium may lead to other effects [1 01 ] . In particu- 
lar, a current drive breaks inversion symmetry, and it is 
well known that symmetry plays a crucial role in criti- 
cal phenomena. In this paper we investigate the effect 
of a current drive on the critical properties of a low di- 
mensional system near a ferromagnetic quantum critical 
point. We present a derivation from microscopies which 
reveals the inversion-symmetry-breaking terms, and we 
analyze their consequences. 

The geometry that will be studied (shown schemati- 
cally in Fig [1]) is a very thin layer laid on top of a sub- 
strate with which the layer can exchange particles as well 
as energy. We take the layer to be a metal which is tuned 
to be near a ferromagnetic-paramagnetic quantum criti- 
cal point. Departures from equilibrium and breaking of 
inversion symmetry are provided by current, j, flowing 
along the layer, while the coupling to the substrate al- 
lows the system to reach a nonequilibrium steady state. 
An earlier paper Q studied a similar system but with 
nonequilibrium provided by current flow across the layer, 
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FIG. 1: A 2d itinerant electron ferromagnet driven out of 
equilibrium by application of an in-plane electric field. A 
steady state is reached via coupling to an underlying sub- 
strate. 



so the issue of inversion symmetry breaking did not arise. 

We now outline the relevant scales and summarize our 
major results. An electron in the layer will escape into 
the substrate after a typical time t sc determined by the 
strength of the coupling to the substrate and taken to be 
long compared to typical electronic times such as the in- 
verse of the Fermi energy Ep\ EpT sc » 1. The escape 
time implies a length l sc = vft sc with vf the Fermi veloc- 
ity in the layer. A current j flowing in the layer implies 
an electric field E = pj. The departure from equilib- 
rium produces decoherence whose consequences may be 
approximately parametrized by an effective temperature 
T e f f [5|, la] which in the present problem we find to be 
given by the voltage drop over an escape length; 



T e ff — eEl s 



(1) 



We present here a theory valid in the limit T e ffT sc << I . 

The new feature of the longitudinal current is a break- 
ing of inversion symmetry. The magnitude of the sym- 
metry breaking effect in the spin sector turns out to be 
related to the spin current ]m = if — j[ induced by an 
applied electric field. In the model we study we find 
that for a magnetization density M the spin current is 
]m = vdM with magnetization drift velocity vd propor- 
tional both to the applied electric field [n], [12] and to 
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deviations from particle-hole symmetry. The importance 
of particle-hole symmetry may be seen in a simple Boltz- 
mann transport picture: the current j a of electrons of 
spin a is j a — J (dk)v(k)n(e k — fx a ;E) with n a distribu- 
tion function depending on the applied electric field, the 
electron energy, and the spin-dependent chemical poten- 
tial. In a particle-hole symmetric situation, v(k) is in- 
dependent of k and £fc = v\k\, so a simple variable shift 
implies that jj = ji implying jm = 0. In order of magni- 
tude vofcp ~ T e ff << Ep, so that the effects of inversion 
symmetry breaking turn out to be small. 

We write a Keldysh path integral 13], 14 1 , make a sad- 
dle point approximation followed by an expansion of the 
Keldysh action in fluctuations about the saddle point [!, 
[l5l | . The resulting nonequilibrium Ginzburg-Landau the- 
ory is then treated by renormalization group [l^ | . The 
Hamiltonian of the system is 



H = Hbath + Hlayer + Hla 



yer—bath 



(2) 



where Hi ayer is the interacting electron layer whose crit- 
ical properties we are interested in, Hbath describes the 
underlying substrate which will be assumed to be an elec- 
tron gas which is always in thermodynamic equilibrium, 
while Hiayer-bath couples the two. Thus, 

Hlayer = ^ C px - ^- c A d l ± a d P± * + U Yl ^ ^ 3 ) 
p ± cr i 

Hlayer-bath = ^ ( tc bath,* d P±<? + h ' C ) ( 4 ) 

with a the spin label, p± the momentum within the layer 
which is a good quantum number, and A = —cEt. The 
only relevant features of Hbath are its density of states, 
Nbath, which implies that t sc = (2TrNb a tht 2 )~ 1 and its 
resistivity, which we take to be very high relative to that 
of the layer so we may couple the electric field only to 
the layer electrons. 

The interaction term is treated via the usual Hubbard- 
Stratonovich transformation, which involves introduc- 
ing Hubbard-Stratonovich fields rrii± in the time-ordered 
(denoted by — ) and anti-time ordered (denoted by +) 
Keldysh axes [5|, choosing a quantization axis at each 
point in space and time 17J and introducing appropriate 
spin rotation matrices R± . To simplify the notation we 
do not write these explicitly except where needed. 

After integrating out the electronic degrees of freedom 
in the bath and layer we obtain the following Keldysh 
action for the classical (m c ; = m -+ m + ) an d quantum 
components (m q — m ~ 2 m+ ) of the auxiliary fields 



Zk = J D \m q m c {D [f2]] i 



— i Jd xdt2Um q m c i 



Trln 



xe 



it GZ L -Ua 



m c i m q 
m q m c i 



(5) 



where the Keldysh Greens function for the layer electrons 



G = ^ G A 



(6) 



The effect of the nonequilibrium drive is expressed via 
the E-field dependence of G, which we now discuss. The 
solution of the Dyson equation for the retarded/advanced 
component of G is 



,-ifti dxe{p+eEx) 



1*1 — *2 I 



G R / A (p,t u t 2 ) = TiO(±(h-t 2 )) 

(7) 

Introducing the time difference t = t\ — t 2 and canonical 
momentum k = p + eET : (T = tl ^* 2 ) we find that up to 
terms of order {eE) 2 {T^ c d 2 ek/ {dk 2 )) in the exponential, 
the retarded/advanced components of the layer Green 
function G take the equilibrium form 

G R/A ( P ,t l7 t 2 ;E) -» G R/A (k 1 r) = TiO(±T)e~ lekT e-^ 

(8) 

By introducing the lattice constant a and noting that 
the fermi energy Ep ~ vf/cl, we may estimate the mag- 
nitude of the neglected terms as 9 e ^p F dk ^e^F^ which 
is much less than unity in the limit of relevance here. 

The Keldysh Green's function G K = G R -f K -fK-G A 
where the distribution function fx = 1 — 2/ obeys 



K 



K 



0*1 



dh 



e p -jL.A( tl )fK + e p -^. A ( t2 )f 



K 



s w • !k + f K ■ ^ 



(9) 



E R _ e a = and S K = _ y, a )(1 - 2g), g be- 

ing the distribution function of the substrate. Fourier 
transforming Eq. [5] with respect to the relative time 
t = *i — *2, changing variables to the canonical momen- 
tum k = p + eET and expanding in E one finds that the 
distribution function at steady state obeys, 



- df 0/ / - de k 
eE ■ -L + [eE ■ 
dk duj V dk 

= — i-f + : g] 

Tan 



i_0V 

24 8lu 3 



eE- 



0_ 

dk 



(10) 



The usual quasiclassical arguments [18| imply that the 
first term in Eq [TU] is negligible while in the the weak 
field limit the third term may be dropped. With these 
simplifications we find / = f s + f a where 



fk,x 

ra 

J k.x 



I x 

siqn(x) T, 5 z * 2 

X) + - V '- e V( eM ^-rsc) ( U ) 



eE ■ v k T s , 



>■{( 



2\ eE ■ v k T s . 



where x — lu — fi and v k = dek/dk. Substitution of Eqs 
[Tl][r2]into EqfTUlthen shows that the neglect of the third 
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term in Eq[lO]is justified at weak layer-substrate coupling 
( E F Jdk3 ^ i E FT sc ) ) while the first term is negligible 
in the weak field limit T eff < E% I 

The quantum critical analysis proceeds [3| by find- 
ing a saddle point with m q = and an m c i that obeys 



2iUm c i 



-UTr 



<TT X 



(13) 



where r are Pauli matrices in Keldysh space, and then 
expanding about the saddle point and using renormaliza- 
tion group methods to treat the resulting action. In the 



disordered phase the saddle point is m c i = 
by expanding the Trln in Eq. [5] we obtain 



m a 



Zk = I T> [m q m c i] e 



-i J d d xdt2Um q m c i 



and 



-(14) 



where 



J K 



-iU 



(HI 



(2n) d 2tt 



[m* q m c iIi R + c.c] + m q m* q [n*Ql5) 



involves polarization bubbles Hr/a/k computed in terms 
of the G. For the Ising case = 0, while = 
J2i=i 4MimJm*p. We have used Eqs [5J [TTIand [T2"lto 
compute the II and ui [5]. Expanding in q, Q, distance 
from criticality 6 and T e //r sc we find 

UIL R (q,n) = 5 + Aq 2 -iT sc (Q-v D ■ q) (16) 

The z(f2 — vd ■ q)T sc factor expresses the non-conserved 
dynamics induced by the escape of electrons from the 
layer into the substrate; the form of this term shows that 
the fluctuations drift at the magnetization drift velocity 
Vd which in the approximations considered here is 



v D = eEr sc —— 
iVo 



d d k 
W) 



1 d(e k - M) 



d 2 e k 
dk 2 



where N — J ^ 2 ^ d S(ek — n)- The drift term was missed 
in previous work [6j . To estimate its order of magnitude 
we note that Eq[17]involves the scale T e f / and bare quan- 
tities. Because one may regard the dimension of velocity 
as [Energy x Length] and the only available length is the 
fermi wave vector, we estimate vd — T e ff/kp. 

The Keldysh polarization bubble contains information 
about decoherence and noise. In contrast to H R / A , H K 
is found to depend on f2 only because the noise arises 
from coupling to the substrate and does not drift in the 
presence of a current in the layer. We find in 2d 



n£(fi) = -2*7 



(18) 



|tl| 



■<> ■- T, rr f ^|cos0|e ^77M 



and for a Id system, 

nf rf (0) = -2 tTs , 



\n\+T eff e T °ff 



(19) 



Eqs [18] and [19] show that the scale T e f f is similar to a 
temperature in that II is proportional to Vl for ft > T e ff 
but is proportional to T e ft for Q < T e ft. An electric 
field induced effective temperature was previously iden- 
tified in Q, where the current flow across the layer con- 
trolled the decoherence, and in Q where the equivalent 
of r sc was argued to arise from electron-magnon scatter- 
ing which itself depended on T e //, leading to a nonlinear 
relation between T e // and E. 

We now present a renormalization group treatment of 
Eq. 1141 The theory involves a momentum cutoff A ~ kp 
and we define the energy unit to be u^A which we assume 
to be >> T e ff. We also assume vp < Vp and u, 5 « 1. 
We integrate out momenta in a small range A to Ae~' and 
simultaneously integrate out all frequencies from — oo to 
oo (this procedure makes the preservation of causality 
manifest [5()- We then rescale momentum q — * qe~ l , the 
frequency u> — -> we _zi , and the fields m, — ► nme al , so as 
to keep the second and third terms in Eq. 1161 invariant, 
implying dynamical exponent 2 = 2. We have 



dT, 



eff 



dl 

dS 

~dl 

dvp 
~df 

dui 

~df 



2T eff (20) 

2S + C lUl (21) 

v D (22) 

[2 - d] Ui + 0(u k ui) (23) 



Here C\ = lim Vf 



C°° dm Un K (w,Te f f) 
u P A/T e// ^oo J_ 00 2tt UIl R (q=A,uj)Un A (q=A,Lj) ■ 

Solving in dimension d = 2 gives T(l) = T° ff e 2 \ 
(17) v D (l) = v a D e l , 5(1) = r e 21 and = 2tt 2 /(9/) af- 

ter initial transients have decayed. Here the superscript 
denotes initial values and ro is related to 8q and u 
as described in [2J. The most interesting regime is the 
quantum critical/renormalized classical one occurring for 

T®jj > |r |. In this case, at the scale It ~ InJl/T^^ a 

crossover to classical scaling occurs. Eq. [22] shows that 
the drift term remains small (ud(Zx)A << T(It) ~ 1) so 



that up to corrections of order ^JT^^/Ep the drift does 
not affect the scaling. The treatment of the crossover to 
the classical regime follows the procedure given in [j| 
leading to a long wavelength theory described by the 
Halpcrin-Hohcnberg model A dynamics, but with an ex- 
tra j ■ V coupling to the order-parameter which represents 
current induced drift. Restoring physical units, 



T,<- [ — VD' V ) //>,./ Ir. / i 



-AV 2 



27r 2 m 2 cl 



9lnJl/T? ff 



m cl + ^ (24) 



4 



where 

(£(x,i)£(x',t')> 



sc 1 eff 
7T 



S(t-t')S(x-x') (25) 



The delta-correlated nature of the noise means that the 
low energy theory is Galilean invariant so that the drift 
may be eliminated by the transformation x ~ * x — vjjt. 

We now turn to the case of a Heisenberg magnet which 
may be treated along very similar lines except that we 
require a vector Hubbard-Stratonovich field 17] which 
we denote by A. The physics is very similar to the Ising 
case except for fluctuations about the ordered state. We 
concentrate on this in what follows. Expanding around 
A c i = (m c i)z, A q — where (m c i) is given by Eq[l3l one 
obtains the Landau-Lifshitz-Gilbert equation of motion 
for the transverse spin wave fluctuations: 



(Al + iBTy) 



with A 



dt v Dxx 



■ V 



(26) 



fluctuations to drift with a velocity related to the spin 
current; however for Ising symmetry we find that in the 
regime where critical fluctuations are important the the- 
ory is effectively Galilean-invariant so the drift may be 
eliminated by an appropriate transformation. In mod- 
els with rotational invariance the current leads to a spin 
wave instability whose critical behavior might be inter- 
esting to study further; however in the models we have 
studied the symmetry breaking effects are numerically 
small so the main effect of the departure from equilibrium 
is an effective temperature, proportional to the voltage 
drop over a distance of order of the mean free path. Ex- 
tension of these results to the superconducting case where 
the direct coupling between the order-parameter and the 
current can give rise to other nonequilibrium effects be- 
sides noise is currently in progress. 
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CV 2 and B = 



-qI — VDxy ■ V — DV 2 where £ is the white noise 
defined in Eq. [53] The spin- wave correlation function ob- 
tained from Eq. [26] may be computed by standard meth- 
ods and is found to evolve in time as e -F '* 1- * 2 ' with 



Re IF] 



(a xx C + a xy D) q 2 + (vDxx — VDxy) 



at 



showing that a dynamical instability signaled by F R < 
occurs if even an infinitesimally small current is ap- 
plied to a magnet of Heisenberg symmetry, provided 
vdxx VDxy For v Dxx = v Dxy , the v D ■ q term may 
be eliminated via a Galilean transformation and the cur- 
rent does not drive a spin-wave instability. vd xx — vjj xy 
is a non-universal quantity. In the model defined by Eqs 
I3l4l f p xx differs from vrixy only in the presence of energy 
or momentum dependent scattering rates. 

A derivation from microscopies of the current induced 
drift within a linear response treatment was presented 
in 1|J 13] for models with impurity scattering. Our treat- 



ment here differs in two ways, one is that we consider 
a clean system where the dominant scattering mecha- 
nism is via an inelastic coupling to an external reservoir. 
Secondly we go beyond linear response, which allows us 
to capture the effect of current induced noise. Ref [21] 
showed that in the absence of noise Eq [^Hl has chaotic 
dynamics in the instability regime; the modifications of 
their results in the presence of noise and critical fluctua- 
tions is an interesting open question. 

In summary we have derived from a fundamental the- 
ory the effect of an in-plane current on a magnetic quan- 
tum critical point. The current flow causes the critical 
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